Random tilings of high symmetry: II. Boundary conditions and numerical studies 
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We perform numerical studies including Monte Carlo simulations of high rotational symmetry 
random tilings. For computational convenience, our tilings obey fixed boundary conditions in regular 
polygons. Such tilings are put in correspondence with algorithms for sorting lists in computer 
science. We obtain statistics on path counting and vertex coordination which compare well with 
predictions of mean-field theory and allow estimation of the configurational entropy, which tends to 
the value 0.568 per vertex in the limit of continuous symmetry. Tilings with phason strain appear 
to share the same entropy as unstrained tilings, as predicted by mean-field theory. We consider the 
thermodynamic limit and argue that the limiting fixed boundary entropy equals the limiting free 
boundary entropy, although these differ for finite rotational symmetry. 

I. INTRODUCTION 

A tiling is a filling, without gaps or overlaps, of a given region of a d-dimensional Euclidean space, with tiles which 
differ according to their shapes, sizes, and orientations. In the present paper, the tiles are d-dimensional rhombohedra, 
which we will generically call "rhombi" in the following. Each tile is the projection of a <i-dimensional face of a D 
dimensional hypercube (D > d) into ci-dimensional space. The difference D — d is known as the codimension of the 
tiling, and we say that we are dealing with D — > d tilings. Details of this construction can be found in a previous 
papesi, hereafter referred to as "paper I" . We mainly focus on the case d — 2. A random rhombus tiling can be viewed 
as a fluctuating membrane in this higher dimensional spaced, the membrane being the union of hypercube faces (see 
paper li, section II. 1). 

Rhombus tilings are dual to de Bruijn gridsi^. In two dimensions, these grid lines pass through the midpoints of 
parallel rhombus edges. Every rhombus edge orientation defines a family of effectively parallel de Bruijn grid lines. 
Dc Bruijn grid lines within a family never cross. In contrast, lines of different families do cross, and their crossing 
defines a rhombus of the tiling. There are D rhombus edge orientations and hence D families of de Bruijn grid lines. 
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FIG. 1: Fixed boundary D — > 2 tilings with D = 3, 5 and 7. The side lengths of the polygons are all p = 8. 

Rhombus tilings provide simple models for quasicrystals^, metal alloys that exhibit rotational symmetry forbidden 
by conventional crystallography. One of the key properties of random tilings^ is their configurational entropy that 
may play a role in stabilizing the quasicrystal state. The source of entropy is localized tile rearrangements known 
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as phason flips. Groups of three adjacent rhombi may be permuted so that their perimeter remains fixed while their 
shared vertex moves to a nearby point. Our present focus, as in paper I, is the limit of high rotational symmetry 
obtained as D — > oo at fixed d. We examine this limit because of its intrinsic interest and because it deepens our 
understanding of finite D tiling models relevant for real alloy systems. 

As usual in statistical mechanics, tiling systems can have different boundary conditions, such as free, fixed or 
periodic. For example, the tilings in figures HI l2~l or PHI have fixed polygonal boundaries. It is believed that free and 
periodic boundary tilings reach equal entropies at the thermodynamic limit. Fixed boundary tilings, in contrast, 
exhibit entropies that are strictly smaller than the free boundary entropy. For example, the fixed hexagonal-boundary 
entropy of 3 — > 2 tilings (see Fig. ^ left) equals 0.261 (— ), while the free boundary entropy equals 0.323 (— ). This 
phenomenon can be understooc&Siifiiii by inspection of Fig. ^ The local entropy density displays a gradient between 
crystalline regions near the boundary, where the entropy density vanishes, and the central region, where the entropy 
density reaches a maximum value equal to the free boundary quasicrystalline value. Only at the very center of the 
tiling does the influence of the boundary disappear. 

Paper li developed a mean-field theory for free boundary tilings, applicable to the limit of high rotational symmetry. 
Two earlier papers^ii presented initial studies of this problem. The first papeu^ proposed an upper bond on the 
entropy in the limit of large D and discussed problems associated with the thermodynamic limit of tiling models. 
The second onei^ presented a preliminary mean-field approach of the entropy calculation. In the present paper we 
focus on fixed boundary tilings because they are easy to simulate numerically. We present results of our simulations, 
including an accurate estimate of the limiting entropy density. We investigate the role of fixed boundary conditions 
and argue for a boundary-condition- independent thermodynamic limit in the limit D — > oo. Similarly, we show that 
phason strain does not influence the entropy for large D. 

The organization of this paper is as follows. We start in section ^ with a discussion of fixed boundary tilings, 
and describe their relationship to interesting problems of pure mathematics. Next, we review the problem of the 
thermodynamic limit in section lTTTl where we argue that free and fixed boundary tilings attain the same thermodynamic 
limit as D — > oo. We also explore finite D corrections to the fixed boundary entropy and their relation to the 
inhomogeneity of the tilings. Then, section llVl describes our Monte Carlo simulations. In that section we explore 
the entropy, path counting statistics, vertex coordination statistics and the role of phason strain. We also confirm 
numerically the results of section ITTT1 The paper is written so that the reader does not need the notations of paper I, 
except in the appendices. 

II. FIXED BOUNDARY TILINGS 

Fixed boundary tilings such as those of figure Q possess the most natural boundary conditions for tilings coded by 
combinatorial objects such as generalized partitions or sorting algorithms. Even though these boundary conditions 
are unnatural for bulk quasicrystals, this inconvenience is counter-balanced by the powerful tools provided by such 
codings in terms of understanding of tiling set structures and enumeration of tilings. Furthermore, we will demonstrate 
in this paper that fixed and free boundary tilings become equivalent in the large codimension limit, which will justify 
a posteriori the use of fixed boundary conditions. 

As already described in the introduction, such fixed boundary tilings lack a proper thermodynamic limit for finite 
D because the boundary has a strong macroscopic effect on the whole tilingLiLliLlLll. This leads to a spectacular 
effect known as the "arctic circle phenomenon-^" in hexagonal (D = 3) tilings, where the tiling is periodic (and 
"frozen" ) outside a perfect circle at the large size limit and random inside this circle. More generally and for larger 
D, fixed boundary tilings have an effective codimension that is smaller near the boundary than in the bulk. Very 
near the boundary the effective codimension vanishes. The tiling becomes a crystalline domain comprising a single 
tile type. The local entropy density vanishes because no phason flips are possible. In the membrane picture, this 
region corresponds to a flat area with large tilt. Further from the boundary the effective codimension grows. The 
variable codimension results in an entropy density gradient, growing from zero at the boundary to the maximal, free 
boundary, entropy density at the tiling center. Consequently, fixed boundary tilings have a smaller total entropy per 
tile than free boundary tilings. 

In contrast, we will show in section ITTT1 that when D becomes large (figure this heterogeneity diminishes and a 
thermodynamic limit is restored. Since fixed boundary tilings are easier to specify and manipulate for arbitrary D, 
they are the appropriate tool to tackle large D entropies or other statistical properties. In the following, the fixed 
boundary entropies per tile will be denoted by ctd and their limiting value by tJoo. Corresponding values without the 
over-bar refer to free boundary tilings. 

In terms of de Bruijn dualization, any two lines of two different families intersect inside a fixed boundary tiling^ 
(the grid is said to be "complete"). Duality associates tile with the intersection of two lines. Therefore if there are fej 



lines in each family Fi , the number of tiles is 



If all h are set to 1, Nt = D(D ~ l)/2. 



N T = k ^ k r 

l<i<j<D 



(11.1) 




FIG. 2: A 40 — > 2 fixed boundary tiling with one de Bruijn line in each family. 

The relationship between tilings and partitions has been widely explored in referencesii^iS^ii. The idea is to code 
a random rhombus tiling by an array of integers satisfying certain ordering constraints. Figure [21 illustrates this 
point in the simple hexagonal casoi there is a one-to-one correspondence between hexagonal tilings filling a centrally 
symmetric hexagon of sides k,l and p on the one hand and sets of integers arranged in a rectangular array k x /, 
decreasing in each row and each column, on the other hand. These integers are bounded between zero and the side 
length p. They decrease monotonically. Such a set of integers is called a plane partition. 
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FIG. 3: Three-dimensional representation of an hexagonal tiling filling a centrally symmetric hexagon (left). This object can 
be seen as a staking of cubes. The height of the different stacks can be arranged in an array (right). They are decreasing in 
rows and columns. This array together with these order relations constitutes a plane partition problem. 

This point of view can be extended to any D — > d problemi^iSiii: any D — > d tiling can be coded with a generalized 
partition on an array related to a D — 1 — > d tiling. This correspondence is also one-to-one. However, tilings coded 
by partitions have fixed polytopal boundary conditions. In two dimensions, the boundaries are centrally symmetric 
2_D-gons (see figure 0. 

We now present an analogy between fixed boundary tilings and algorithms for sorting lists**. Although Computer 
Science is not our motivation, we are interested in the same enumeration problems***—. More generally, the rich 
topological and combinatorial properties of random tilings make them an active field of research in pure mathematics, 
combinatorics and computer sciencei2i22i2ii22. We shall use results from these fields throughout this paper. 

In the sorting language, a comparator [i; j] acts on a list (xi, x%, ■ ■ ■ , Xb) of numbers as follows: X\ and Xj are 
respectively replaced by min(a;j, Xj) and raax(a;,, Xj). Following Knuthi&, we call a complete sorting algorithm a 
sequence of such comparators which sorts in the increasing order any list of real numbers (x\,X2, ...,Xr>)- This 
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sorting algorithm will be called primitive if each comparator can be written [i;i + 1]. We also suppose that this 
algorithm is not redundant, that is to say it does not contain any comparator that could be suppressed because 
previous comparators already insure that Xi < Xj . Knuth shows that a sequence of comparators is a sorting algorithm 
if it correctly sorts the completely reversed list [D,D — This means that a complete primitive sorting 

algorithm is a sequence of comparators [i; i + 1] that transforms the list (D, D — 1, . . . , 1) into the list (1,2,..., D). 

Such an algorithm has a diagrammatical representation in which the D variables Xi are represented by D horizontal 
lines. Each comparator [i;i + 1] is represented by a crossing between lines i and i + 1. Figure 2] illustrates this 
construction. A continuous line follows a number during the sorting process. Since every number must be compared 
to every other number, and since there is not any redundancy, every line crosses every other line, only once. There 

are ( ^ J crossings. 




FIG. 4: A diagram associated with a sorting algorithm acting on five element lists. A line follows a number during the sorting 
process. Each pair of lines cross, only once. Each crossing represents a comparator [i; i + 1]. 

We now establish the link between sorting algorithms and the de Bruijn representations^ of D — > 2 tilings. Each 
line of the diagram ^represents a de Bruijn line and crosses every other line exactly once. However, different sorting 
algorithms sometimes represent the same de Bruijn grid, since only the crossing topology is meaningful. For example, 
in figure the fourth and the fifth comparator (i.e. [4;5] and [2;3]) are applied in this order (these comparator are 
circled in the figure). If they were applied in the reverse order, the algorithm would be different whereas the de Bruijn 
grid would be the same. 

Therefore we define equivalence classes of sorting algorithms 18-21 . We say that two successive comparators [i; i + 1] 
and [j; j + 1] commute if \i — j\ > 1. Two algorithms are equivalent if they differ by a finite number of comparator 
commutations. Equivalence classes of D-element sorting algorithms are in one-to-one correspondence with Z?-family 
grids with one line per family, and therefore with tilings inscribed in polygons of side 1. Following Knuth, we denote 
by Bo this number of equivalence classes, whereas the number of algorithms is denoted by Ar> {Ar> > Bjj). Since 
each pattern of crossings (equivalence class of algorithms) defines a tiling, the tiling entropy density per tile is 

a D = —logB D . (II.2) 
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FIG. 5: Entropy per tile od of fixed boundary tilings (p — 1). We indicate the limiting values 0.568 (full line) and the close 
upper bound log 2 ~ 0.693 (dashed line). 

Sorting algorithms are easy to manipulate and provide efficient enumeration numerical tools. This analogy provides 
us with the number Be, of D — > 2 tilings for small values of D < 10 from the work of Knuthi& as extended by usi£ 
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(see Table IJ). Unfortunately it is not possible to reach the next value (D — 11) for the foreseeable future using our 
current algorithms. Figure [S] suggests that the entropy (Joa{p — 1) has a well defined limit when D tends to infinity. 
Indeed, Knutb»& gives lower and upper bounds on Bp from which we get 

ilog2< lim <7£,<21og2. (II.3) 

3 D—>OG 

Moreover, Bjorner (see reference^, p. 270) derives a better upper bound: 

lim 5 D < 1.44 log 2. (II.4) 

In the following section, we argue that large D, fixed boundary tilings, have the same entropy per tile as free boundary 
ones. Since we demonstrated in paper I* that for free boundary entropies lim^^oo od < log 2, we finally get a better 
bound: 

lim a D < log 2. (II.5) 

This point of view generalizes to systems with more than one line per de Bruijn family 15 , leading to the definition 
of partial sorting algorithms for merging pre-ordered lists of numbers. Suppose that we have D families of ki numbers 
each (i = 1, . . . , D), and that in each family the numbers are presorted in increasing order. Algorithms which order 
the union of these sets of numbers are called partial sorting algorithms. The ideas are essentially the same as in the 
previous case, except that, since the numbers of a given family are already ordered, the corresponding lines do not 
cross within the family. The corresponding diagram is similar to a de Bruijn grid with D families of lines, ki lines in 
each family The tilings are equivalence classes of such algorithms. They are inscribed in polygons of sides k\, . . . , hp. 
In reference (— ), we derive an analytic expression for Ao(ki, . . . , ko)- However, it proves impossible to calculate Bd 
analytically. 

The diagonal case, with k\ = fca = • • • = fco = p will be of special interest below. Let Bd(j>) denote the number of 
such tilings. Using the sorting algorithm analogy, we computed the numbers of tilings inscribed in polygons of side 2 
(B D (2)) for small values of D (see Tabled. 



III. THERMODYNAMIC LIMIT 



In this section, we argue that the fixed boundary entropy in the large D limit equals the free boundary entropy and 
is independent of the phason strain. Proofs are provided in appendices ^ and [5] Free boundary and fixed boundary 
entropies can be compared when both are known. For hexagonal tilings with D = 3, the free and fixed entropies are 
respectively 0.323& and 0.261—. For octagonal tilings with D — 4 the values are 0.434^ and 0.36(l)i 5 ". Therefore the 
relative difference (a — a) /a decreases between D = 3 and D = 4, consistent with a vanishing difference at large D. 
Exact data on p = 1 and p = 2 tilings (Tables [I] and [H} suggest that when D becomes large, the entropy becomes 
independent of p. 

A glance at figure ^ suggests that when D becomes large, the faceted regions of fixed boundary tilings, which are 
caused by the strong influence of the boundary^, occupy a vanishing fraction of the tiling area close to the boundary. 
In these regions, few families of de Bruijn lines cross, so only a small number of tile types are present. The central 
region of the tiling, where a large number of de Bruijn line families cross, becomes dominant as D increases. The 
coarse-grained entropy density in the central region approaches the large D free boundary entropy a^, . 

More precisely, let (xi) be any set of numbers < a < Xi < (3, where a and (3 are fixed positive real numbers. 
Consider fixed boundary tilings of a 2D-gon with edge lengths ki = Xip as p — > oo. We impose such bounds on the 
side lengths in order to be sure that our tilings are effectively large codimension ones (see also paper li (section II.3)). 
The tile fractions are 



lJ N T kik 2 + .., + kD^kD 1 

owing to eq. (|II.1|I . In an unstrained tiling, where all ki are equal, 

2 



(III.1) 



« D{D - 1) 

These quantities differ noticeably from their free boundary counterparts^. The condition < a < Xi < (3 imposes 

2 

/3 2 n ' l] ~ P 2 D(D - 1 



(III.2) 



Uij > n* 3 = — — . (III.3) 
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This lower bound ensures that no tile fraction is vanishingly small and therefore that the tiling is really a large 
codimcnsion one rather than a small codimension one with few defect lines added in. 

Under this condition, appendix E] proves that the sequence od reaches a thermodynamic limit Ooo that does not 
depend on the set (xi). Furthermore, these tilings are locally equivalent to free boundary high symmetry tilings. The 
local entropy density thus approaches <7qo, and Oqo = lirnr>— >oc <?d — &oo- Moreover, this result remains valid if D 
becomes large while p is held fixed (even in the extreme case foj = p = 1). In particular, the limits D — > oo and p — ► oo 
commute with each other. The proofs of all these results are given in appendix ^ This insensitivity of entropy to 
boundaries could be anticipated since we know that phason elastic constants vanish at large D fi] and that, because 
of the regular character of large D tilings established in appendix [5] the phason gradient E is nearly everywhere 
bounded. 

We find in appendix^] that large D tilings have a central region holding half of the tiles, which contains lines of all 
de Bruijn families and all tile species and which is nearly homogeneous and strain-free. Away from this central region, 
the tiling no longer contains lines of all de Bruijn families. We define the effective (coarse-grained) dimension D e ff 
as the number of de Bruijn families present in a small tiling patch, and the effective codimension c e // = Deff — 2. 
This effective codimension decreases from D — 2 at the center to at the boundary. For diagonal tilings (xi = 1 so 
that ki = p for all i), where the tiling has a circular symmetry at the large D limit, the effective codimension c e // 
varies as 

c e //(r) = (£>- 2) 7 (r/R) (III.4) 

where r is the distance to the center and R is the radius of the tiling. The function "/(r/R) varies from the value 
7(0) = 1 at the tiling center to 7(1) = at the boundary. For large D, c e //(r) diverges for all r < R, so the local 
entropy density approaches Oqo for all r < R. 

This dependence of c e / / on the location in the tiling can be used to estimate finite size corrections to or> ■ Consider 
a fixed value of D 3> 2, and find the number of tiles in the regions where c e // < Co, with Co <C D — 2 some fixed value. 
These outer regions form an annulus R an < r < R. Inside the radius R an , the mean entropy density is nearly doo, 
but in the annulus the mean entropy density is o an < ojj. We can estimate the mean entropy density on the whole 
tiling of radius R as 

ob ~ no an + (1 - n)(Too, (III.5) 

where n is the fraction of tiles in the annulus. 

Now it is shown in appendix equation (|A.16Jl , that the fraction of tiles in the region of effective codimension 
c e ff = xD is n(x) — ip(x)/D where ip is a regular function tending to 2/3 when x goes to 0. Hence the fraction of 
tiles in the annulus 

-55- < n «> 



Finally, we estimate 



B 

D ~ (Too - — , (111.7) 



where B > is a constant related to the entropy difference (ooo — o an ) and the chosen value co, both of which are 
independent of D. 

This result holds only for a large tiling side length p because the annulus width must be large as compared to a tile 
edge in order to define properly an entropy per tile. Appendix IA 31 savs that the effective codimension c e // = cq at 
Ran if Ran = (1 — e)R where 

^2 „1 

(HI.8) 



24 (D - 2) 2 

The above requirement reads eR 3> 1, that is to say p 3> D with R = pD fir . The large p limit must be taken before 
the large D one. Numerical finite D corrections in small p tilings turn out to be of order 1/D but positive (see next 
section). 

In conclusion, the local structure, and therefore the entropy per tile, of large D tilings is independent of their size, 
shape, tile fractions and boundary conditions. In all these cases, the local structure of these tilings is similar to free 
boundary large D tilings, where one encounters only one de Bruijn line per family in large tiling patches^. In this 
sense, it is a true thermodynamic limit. These results will be confirmed by numerical simulations in IIV 31 

Henceforth, we concentrate on the study of tilings inscribed in 2D-gons of sides ki = p = 1 as representative 
of the whole class of high symmetry tilings. In particular, we focus on p — 1 tilings to extract (numerically) the 
thermodynamic limit of the entropy density. 
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IV. SIMULATIONS 



Monte Carlo numerical simulations^* are widespread in the fields of random tilings and quasicrystals where there 
exist many difficult unsolved theoretical questions^*. To check our claims of a thermodynamic limit, to obtain Bp for 
D > 10, and to get more precise numerical values than those obtained via mean- field arguments*, we perform Monte 
Carlo calculations on large D tilings. The D — > 2 configuration space is sampled for large D systems via single-vertex 
flip dynamics. This method is validated by the connectivity of the space of configurations in two dimcnsions*£i22*2i. 

These simulations utilize fixed boundary tilings for the sake of technical convenience. It is easier to code a fixed 
boundary tiling in the memory of a computer, and the entropy is easier to estimate in this case. We have established 
that the central region of a large D tiling is close to a free boundary tiling. All the statistics related to free boundary 
tilings will be collected in such central regions. Note that, as far as unstrained tilings are concerned, we restrict our 
study to polygonal boundaries of side p — 1 since we establish that such tilings behave like large p ones when D goes 
to infinity. 



The algorithm is implemented as follows: at each Monte Carlo step, select a vertex at random, with uniform 
probability; if this vertex is flippable, then flip it; repeat sufficiently many times. The method is validated by the 
ergodicity of the space of configurations via flip dynamics: any configuration is reachable from any initial configuration. 

A key point is to check that this algorithm samples configurations with uniform probability, since we are only 
interested in configurational entropy in which all configurations play the same role. Thus we need to establish that 
this algorithm defines a Markovian process which satisfies the detailed balance conditior 24 : 



where w(C\ — > C2) denotes the transition probability from configuration C\ to configuration Ci. When two configura- 
tions differ by a single flip, there is only one way to go from one of them to the other one by a single flip. Therefore 
W = l/Nv where Ny is the number of vertices. Since Ny is independent of the tiling, relation (|IV.ljl is established. 

Note that the flip acceptance rate equals the fraction (about 31% as reported in section II V 3 1|> of vertices with 
threefold coordination. Vertices with more than threefold coordination cannot be flipped by an "elementary" phason 
flip. An alternative Monte Carlo dynamics, in which a threefold vertex is flipped in every step, does not obey detailed 
balance because the number of flippable vertices is not conserved, and hence eq. H1V.1|) is violated. 

Another point to consider in Monte Carlo simulations is the correlation time to (in units of Monte Carlo steps 
per vertex) which measures how many steps are necessary between samples to avoid excessive sample-to-sample 
correlations. Even though some recent work brought new results for small D (see2* and references therein), this 
question has no rigorous definitive answer for large D. However, our numerical estimations of to are all in agreement 
with the conjecture to < iVy/2. 



To estimate the entropy <jd of fixed boundary tilings, we shall use a path-counting algorithm based upon de Bruijn 
line enumeration^ directly derived from the general method described in section 2.2 of paper li. Figure illustrates 
the method. At each step, we build a D + 1 — > 2 tiling from aD^2 one. Starting from any D — > 2 tiling, we 
choose a path of length D along tile edges, going from bottom to top. We cut the tiling along this path, separating 
the two parts by length 1, then draw new bonds connecting previously identical vertices. Finally, we adjust all edge 
orientations to match the set of symmetry D + 1 . 



1. Monte Carlo algorithm 



w(Ci — ► C 2 ) = w(C 2 — » Ci), 



(IV.l) 



2. Path counting and entropy estimate 




FIG. 6: Iterative construction of fixed boundary tilings (D — ► 2, p = 1). 
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Let Pd denote the mean number of bottom-to-top paths on the D — > 2 tilings: 

p d = w J2 Pd ^ ( IV - 2 ) 

T— 1 

where Pd(t) is the exact number of paths on the tiling r, and Bo is the number of such tilings. Then 

B D +i = PdB d . (IV.3) 
Iterating relation (|IV.3|) , and taking the logarithm yields 

D-l 

\ogB D =J2^gP k . (IV.4) 

k=2 

For large D, the ratio jj log Pd approaches a finite limiting value. Taking the limit of i|lV.4[) as D — > oo, and noting 
eqs. (III. Ill and 1)11. 2|) . we find 

lim an = hm — logPc- (IV. 5) 

Hence, by accumulating statistics on the number of paths allowed on tilings, we may evaluate the entropy. 

We use Monte Carlo sampling to generate an ensemble of D — > 2 tilings. On each tiling r we can quickly evaluate 
-Pd(t) using a generalization of the Pascal Triangle construction: starting at the bottom of the tiling, assign each 
vertex an integer value equal to the number of paths that reach it from the bottom. The value at any vertex is 
iteratively the sum of the values at each prior vertex to which it is connected (see Fig.[7J|. At the end of the process, 
the value of the top vertex is Pd (t) . 




FIG. 7: Path count construction, as described in the text. 



We tested this algorithm on small D systems, where numerical values are in good agreement with the exactly 
known mean number of paths Pd = Bd+i/Bd (table HJ. For example, table ITTT1 shows the convergence towards the 
exact value for D = 9 tilings. For large D systems, simulations are our only means of obtaining information about 
path count statistics, which are displayed in figure|Hl In this figure, error bars are calculated as follows: the standard 
deviation A is measured numerically. The error is estimated by A/y/l, where I is the number of independent measures 
in the simulation (that is^ its total length divided by twice the correlation time to)- Note that the actual errors 
calculated in table ITTT1 are compatible with the previous error bar. 

This graph presents a maximum near D = 20 which is related to finite D and finite size corrections. Paper I 1 shows 
that for free boundary systems, finite D corrections to log Pd/D are of order 1/D 2 . Section ITTT1 of the present paper 
argues for negative 1/D finite size corrections with fixed boundary conditions in the large p limit. Here, though, we 
examine p = 1 systems with one de Bruijn line per family. Our numerical results indicate that corrections also fall 
off as 1/D to the first power. Combining these corrections suggests 

The numerical results in figure [H] reproduce this behavior. All numerical data beyond D — 15 coincide with the 
fitted ones up to error bars. In particular, the maximum observed near D = 20 is reproduced. We estimate the 



9 




FIG. 8: Value of log Pd/D as a function of D up to D = 120; inset: large D behavior as a function of 1/D and quadratic fit. 
The numerical values (symbols) as well as the fitted ones (lines) present a maximum near D = 20. 



limiting value of log Pd/D = 0.5676 ± 0.0001. The uncertainty estimate on this last limit comes from excluding the 
data with D = 120. These observations provide the first order correction to the entropy at p = 1: relation l|IV.4|) 
reads 



D-l 



log Br, (1) ~ fc ( CT ~ + T) - D{D 1) fT °° + A(D - 2 ) ( IVJ ) 

and 



fc=2 



2 2A 
a D (p = 1) = — - _ log^(l) ~ a x + — . (IV.8) 

Note that our simulations yield A > which means that these corrections of order 1/D are positive. By contrast, 
corrections to in relation l|III.7|) are also of order 1 / D but negative. The sign of these large p corrections is the 
combination of two effects: the paths visit low codimension regions near the boundary which should increase the 
entropy (for example, a path running on a square grid of codimension has an entropy density of log 2 > Uoo)\ the 
entropic repulsion between paths near the top and bottom vertices, where they are crowded, decreases the entropy. 
In the p = 1 case, paths are not subject to entropic repulsion. Presumably the coefficient A in eq. 1IV.6(I is a function 
of p and changes its sign when p grows. 



3. Thermodynamic limit revisited 

Here we test the results of section IIIII where it was demonstrated that the local structure and the entropy per 
tile of tilings becomes independent of boundary conditions when D becomes large. We compare relevant numerical 
quantities both in the whole tiling and in the central region, that is supposed to be close to a free boundary tiling. We 
concentrate on vertex coordination numbers which are indicators of the local structure and on path counting which 
is related to entropy. In the last subsection, we focus on strained fixed boundary tilings. 



1. Vertex statistics 



Other pertinent information available from Monte Carlo simulations concerns vertex statistics. They are indicators 
of the tiling's local "microscopic" structure. They also are a fundamental ingredient of the mean-field approach of 
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paper li. We focus on quantities which can be computed in the mean-field theory, especially coordination numbers 
and related quantities. We gather statistics both in the whole tiling (excluding vertices actually on the boundaries) 
and in a central region containing about 20% of the vertices. Presumably, the central region is only weakly influenced 
by the fixed boundaries, and we take that data as representative of free boundary tilings. 



Z=3 



Z=5 
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0.315 




0.31 
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0.01 0.015 
1/D 



0.025 



FIG. 9: Fractions of Z = 3 and Z — 5 vertices in D — » 2 vertices as a function of 1/D, and quadratic fits. Diamonds concern 
the central region (20% vertices), whereas circles concern the whole tiling. The extrapolated large D values coincide. 



We have run Monte Carlo simulations up to D = 100. After plotting the data as functions of 1/D, we extrapolate 
the limiting values via quadratic fits. Figure [§] illustrates two examples of coordination number statistics. We find 
the fractions of vertices (in the whole tiling and in the central region) whose coordination numbers Z = 3 and the 
fractions whose coordination numbers Z — 5. For finite D we see the central region has a relatively small fraction 
of tiles with Z = 3 in comparison with the whole tiling, and a relatively large fraction with Z = 5. The boundary 
regions are thus more likely to have Z = 3 and less likely to have Z = 5 than a free boundary tiling. However, the 
extrapolated values for D — > oo agree to within the accuracy of the extrapolation. The same conclusion holds for 
Z = 4 (not represented here). For larger values of Z, the fractions of vertices are too small and we did not obtain 
relevant measures. Recall that the mean value of Z is exactly 4 in an infinite tiling, according to Euler's theorem. 

Hence, in the limit D — > oo the whole tiling exhibits coordination statistics similar to a free boundary tiling. This 
point supports the existence of a thermodynamic limit - the local structure becomes uniform throughout the tiling, 
independent of proximity to the boundaries and the local strain they create. 

Other quantities we examined relate directly to mean-field theory and were reported in paper li. Our mean-field 
theory was based upon the number of choices to be made while inscribing directed paths on tilings. The number of 
choices N c (v) at a vertex v is the number p of "arms" emerging from a vertex. The number q of "legs" leading into a 
vertex together with the number of arms obey p + q = Z. Paper I reports the simulated probability distribution for p 
and mean value of pq, and shows these agree well with mean- field theory. Here, we display in table llVl the probability 
distribution for Z, and compare it with the predictions of mean-field theory. The agreement is also satisfactory. 

An intriguing feature of the fixed boundary tilings is the divergence of tile vertex density near the boundaries, 



d(r/R) 



1 



a/6(1 - r/R) 1 



(IV.9) 



caused by the vanishing of tile area, and established in appendix IA 31 for a particular tiling. This effect is plainly 
visible in Fig. ^ Given this spatial nonuniformity it is natural to worry about the path counting arguments because 
the choice statistics might vary among different portions of the whole tiling. It turns out this does not happen. Fig. 1101 
illustrates this point by plotting vertex statistics averaged over different portions of the tiling. We divided tilings into 
20 concentric circular bins and evaluated in each bin the average vertex density and the average fractions of vertex 
with each coordination number Z. The results are illustrated for the cases D — 60, 92, 120. Clearly the vertex density 
converges to a non-constant function which diverges at the tiling boundary according to the prediction of eq. I|IV.9|) . 
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However, the fractions of each vertex type rapidly converge to their large D limit independent of position within the 
tiling. Consequently the entropy per vertex should be uniform, even though the entropy per area diverges due to the 
diverging vertex density. 



Radial Variation 




0.2 0.4 0.6 0.8 1 



FIG. 10: Radial variation of vertex statistics from tiling center to edge. Only the vertex density shows spatial variation. Inset: 
vertex density for D=200 (points) matches the function 1/ ^/6(1 — r/R) (line) as r/R — > 1. 



2. Path counting in the central region 

To confirm that the entropy per tile (or vertex) is independent of boundary conditions, we repeat calculations like 
section llV 21 but we concentrate on the central region. It is reasonable to suppose that statistics in the central region 
should match the statistics inside free boundary tilings. 




FIG. 11: Paths in the central region going from a vertex in layer D/A close to the vertical diameter to any vertex in layer 3D/4. 



On a D — ► 2 tiling we define D + 1 layers of vertices as follows: layer contains only the bottommost vertex u , 
layer 1 contains the vertices at distance 1 from vq , layer k contains the vertices at distance k from vo when one follows 
bottom-to-top paths. Layer D contains only the topmost vertex vd- Let P' d denote the mean number of length D/2 
paths from a given vertex v sta rt of layer D/A to any vertex of layer 3D/4. For convenience, we take D a multiple of 
4, and also require that v s tart is close to the vertical diameter (see figure ITTjl . Most paths stay close to this central 
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diameter because their typical deviation from this diameter grows like v D whereas the size of the tiling grows like 
D. Consequently, P'd counts paths in the central region. Data is averaged over many independent tilings. 
Fitting finite D data as in section llV 21 we find the limiting value 

lim l °^f' D = 0.5670 ± 0.0005. (IV. 10) 

D— too D/2 

This value coincides with its fixed boundary counterpart up to error bars. If the path counting statistics in large D 
fixed boundary tilings were influenced by boundary conditions, we would have expected the central value to differ 
from the value over the entire tiling. We conclude that the fixed boundary entropy per tile matches the free boundary 
entropy in the D — > oo limit. 



3. Phason strained tilings 

The data just reported in section llV 3 II suggests that the tiling structure does not depend on the phason strain, 
since statistics are identical in the central region and in the boundary vicinity. We want to test this point directly by 
performing simulations on phason strained fixed boundary tilings. 

We consider three different forms of phason strain in our numerical simulations: low frequency strain in which D/2 
consecutive line families occur with fcj = 1 (i = 1, . . . , D/2), the remaining D/2 consecutive line families occur with 
ki = 2 (i = D/2 + 1, . . . , D); the same low frequency with a larger amplitude so that ki = 3 (i = D/2 + 1, . . . ,D) 
replaces ki — 2; high frequency strain with fcj = 1 (i odd) and ki = 2 (i even). Figure IT21 illustrates an equilibrated 
tiling with applied phason strain. The tiling has D = 30 and the strain is low frequency and large amplitude. 

Computer simulations of path count statistics for D = 10, 50 reveal no significant dependence of path count on 
phason strain amplitude or frequency. The small variations seen (see table EJ cannot be separated from finite size 
effects. Thus entropy does not depend on phason strain, supporting a posteriori arguments developed for free boundary 
tilings 1 . Recall that strain-independence of entropy in fixed boundary tilings relies upon strain-independence in free 
boundary tilings (see section ITTT|) . 



FIG. 12: Phason strained large D tiling. The added paths still go from bottom to top. 

The same conclusion holds concerning vertex statistics on coordination number and choice distributions. No signif- 
icant differences are found in comparison to unstrained fixed boundary tilings or free boundary tilings (see figure IT^)) . 
These data support the existence of a universal local structure, independent of strain and boundary conditions in the 
limit D — > oo. 



V. CONCLUSION 



This paper tackles random tilings of high symmetry with fixed boundary conditions. In random tiling theory, 
boundary conditions are crucial because finite codimension fixed boundary tilings have a lower entropy than free 
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FIG. 13: Fractions of Z — 3 and Z = 5 vertices in strained tilings (p — 1, 1, . . . , 3, 3, . . .) (triangles), together with previous 
data on fixed boundary unstrained tilings (circles) and free boundary ones (diamonds). The limiting values coincide with the 
previous one, up to statistical errors. Lines are quadratic fits. 



boundary ones. However, we argue that boundaries become irrelevant when the codimension becomes large. We also 
demonstrate that p = 1 tilings with one de Bruijn line per family have the same entropy per tile as tilings with p 
large and even infinite, for sufficiently large D. 

As a consequence, the numerical study of the entropy of large codimension random tilings can be concentrated on 
fixed boundary tilings filling a 2D-gon of side lengths set to 1. We perform exact enumeration for D < 10, thanks to 
an analogy between fixed boundary tilings and some class of sorting algorithms. We use Monte Carlo simulations for 
larger tilings, up to D = 120. In both cases, the fact that tilings have a fixed polygonal boundary greatly simplifies 
their encoding and their manipulation in the memory of a computer. Our Monte Carlo data analysis is based on 
the same iterative process (i.e. concentrating on values of Pd) as our mean-field theory in paper I 1 . We obtain a 
very accurate estimate of the entropy per tile, CTqo = 0.5676 ± 0.0001, compared to the approximate mean-field value 
a^ F ~ 0.598. 

Insensitivity of the entropy per tile and vertex statistics to boundary conditions at large D restores a thermodynamic 
limit which does not exist at finite D. However, caution is required. Even if "topological" quantities such as entropy 
and vertex statistics become homogeneous in the tiling, its "metric" properties remain heterogeneous, as discussed 
in section IIV 3 II Consequently we find the entropy per vertex is homogeneous while the entropy density per area 
diverges near the tiling boundary. 

Many of our methods and results can be generalized to higher dimensional tilings. However, the generalization may 
not be simple to implement in practice. Indeed the de Bruijn directed paths become directed surfaces built of faces 
of rhombohedra. The description of these surfaces in terms of successive choices must be generalized, which highly 
complicates their enumeration. 

The restoration of the thermodynamic limit in the high symmetry limit should remain valid in three- or higher- 
dimensions, because the vision of large regions with at most one de Bruijn surface per family still holds. The 
explicit proof should be similar to the appendix [S] with only minor changes to account for the possible spontaneous 
decomposition into regions of lower effective codimension (as in the "arctic octahedron", see RefJ4,). Presumably, the 
effective codimension still goes to infinity for almost all tiles. 

One of the main motivations for the study of high codimension random tilings was the hope that mean-field theory 
might become exact in the high codimension limit. The principal conclusion of this numerical study together with 
the mean-field results of paper I, is that the mean-field theory does not become exact, at least in the simplest version 
which neglects vertex correlations. Rather, the high codimension limit is quite nontrivial, and an exact analytic 
solution is yet to be found. 
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APPENDIX A: THERMODYNAMIC LIMIT OF TWO-DIMENSIONAL FIXED BOUNDARY TILINGS 



In this technical appendix, we demonstrate that the fixed boundary entropy dr> attains a finite limit Uoo when 
D — > oo, and that this limit coincides with the free boundary entropy (Toq. Moreover, we argue that this limit is 
shared with finite side polygonal boundaries as well as non-diagonal tilings. 

Our demonstration relies on the variational principle introduced in references [IBill] to characterize the entropy 
gradient between the center and the boundary in typical tilings with given fixed boundaries. Typical tilings are those 
which maximize an entropy functional defined as the integral over the tiling of a local entropy per tile; This local 
entropy is the /ree-boundary entropy calculated with the local fractions of tiles. The entropy per tile of fixed boundary 
tilings is the maximum of this functional. 

We show below that, given any domain in the tiling, when D goes to infinity, the tiling in this domain is a piece of 
(free) D' — * 2 tiling, such that D' goes to infinity as D does. Hence the local entropy per tile tends to the large D free 
boundary entropy (Too (nearly) everywhere and fixed boundary tilings have the same entropy as free boundary ones. 

We first consider the case where the sides of the polygonal boundary share the same length p (diagonal tilings), and 
take the large p limit before taking the large D limit. Secondly, we consider the case were D becomes large at fixed 
p. We also discuss the non-diagonal case where these side lengths might be different. Our presentation is heuristic 
but can be made rigorous. Finally, we examine scaling laws for several quantities studied in the paper, such as the 
effective codimension, when they are written as functions of the radius r from the tiling center. 



1. Diagonal tilings 

Our demonstration uses a particular tiling denoted as T . It is the dual of a particular de Bruijn gridi^, consisting 
of D families of p straight lines equally spaced out and forming a regular fanSI (see figure EI left). The de Bruijn 
families are denoted by F\, F%, ■ ■ ■ , Fp, and are labeled counterclockwise so Fk makes an angle kit / D with an arbitrary 
reference direction. Figure El displays such a grid and its dual tiling. For finite D, such a tiling is knownAi not to 
maximize the entropy functional among the tilings with the same boundary (see below). By construction, the tiling 
To is divided into domains where only a fraction of the de Bruijn families intersect. In each domain, the tiling is 
homogeneous, with a well defined local entropy per tile at large p. Since only D' < D de Bruijn families intersect in 
such a domain, this entropy a priori differs from the free boundary D — > 2 one. However, by rotational symmetry 
all the domains where exactly D 1 families intersect have the same local entropy. In the following, the union of such 
domains where only D' families meet will be denoted by Ad 1 and called the region of "effective dimensionality" D'. 

Define S®^ 2 as the entropy functional evaluated at To. It satisfies Sq~* 2 < <jd- Then 

lim S^ 2 < a x < aoo. (A.l) 

D— >oo 

Next we intend to show that for any real number h < 1, 

lim S^ 2 >haoo. (A.2) 

D— >oo 



It then follows that 



and, owing to relation (|A.1|I . that 



lim S^ 2 > (A.3) 



aoo < lim S,f~* 2 < CToo < CToo, (A. 4) 

D— >oo 



and therefore that the 3 involved quantities are equal. 

Let us prove the above statement (|A.2|) . We use the fact (see section HV 3 31 and also paper li), that the entropy 
of free boundary tilings does not depend on strain in the large D limit. We show that the local entropy of Tq equals 
CToo in (nearly) all regions Ad 1 ■ 



15 




FIG. 14: The de Bruijn straight lines (left) and the corresponding rhombus tiling in the 4 — > 2 case with p = 20 (for clarity's 
sake, we have only drawn 10 lines per family instead of 20 on the left figure). The regions where only 3 or even 2 families 
intersect clearly appear on the left figure. The corresponding regions are delimited by the two internal octagons on the right 
figure: in the outer region A2, only 2 families meet; In the intermediate one ^3, 3 families meet; And in the central region Aa, 
all 4 families intersect. 

First we estimate how many tiles each region contains. Since a tile is defined as the intersection of two de Bruijn 
lines, we will count the number of such intersections in a given region, by calculating the number of intersections per 
unit area in each region on the one hand and the area of the regions on the other hand. However we work not in the 
"tiling metric" but in the "grid metric" instead. That is to say, in the figure where the de Bruijn lines were originally 
drawn straight (left-hand representation in fig. ll4|) . Figure^l(l e ft) displays a schematic representation of a grid were 
the different regions under consideration clearly appear. 

The number of intersections in a region A of area a is the sum over all pairs of de Bruijn families of the number 
of intersections of a pair. Consider a given pair of families Fi and Fj, i < j, which make an angle 0y. The number 
of such intersections is the area a divided by the area ao of the unit cell of the lattice defined by these (only) two 
families (figure ITSI right). If the distance which separates two lines of a family is set to 1, then ao = 1/1 sialyl- Hence 
the total number of intersections is 

where the indices i and j run over the families present in A. 




FIG. 15: Left: Schematical representations of the D families of de Bruijn lines and of their intersections, to be compared to 
fig. HH , l e ft- When D tends to infinity, the polygonal boundary tends toward a circle of circumference 2Dp and the polygonal 
central region Ad tends toward a circle a diameter p. Right: Two families of de Bruijn lines making an angle 9ij. The number 
of intersections of these two families in a region equals its area a divided by the area ao of the grayed unit cell. 



As displayed in figure [TBI (left) a region Ad 1 is a crown made of 2D kites (or triangles for the inner crown) of 
equal areas. A kite (or triangle) of Ad' will be denoted by K.D' and its area by kd>. Henceforth the area of Ad 1 
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is equal to old 1 — 2Dkd'- In !Cd>, there are D' adjacent families of lines, for example F±,..., Fq> . The angles that 
they make relatively to an arbitrary reference direction are therefore (9; = 9q + It:/ D, I = 1, . . . , D' and the number of 
intersections in JCd> is 

Qd> = Kd,\ >: sin^-^l. (A.6) 

When D goes to infinity, x = D'/D becomes a continuous variable which represents a "fraction" of the D families 
of lines. For example, the previous sum in eq. I|A.6|) can be estimated by an integral on the variables y = j j ' D and 
z = i/D. If E(x) denotes the integral part of x, then 

(j — i)n f x f v D 2 

4>d{x) = 2_^i snl Jj — D 2 j dy j dzsmn(y — z) = — j [nx — sin %x\ . (A. 7) 

l<i<j<E(xD) 



JO 



The central region Ad plays a particular role among the regions Ad 1 '- when D is large, the fraction of tiles in 
Ad tends to 1/2. Indeed, as illustrated in figures ITU Cleft) and 1151 Cleft), the central region tends toward a circle of 
diameter p. Its area is irp 2 /4, and, owing to the above result applied to the case x = 1, the number of intersections 
that it contains grows like 

D 2 irp 2 D 2 p 2 

- ~r ■ (A - 8) 

D 2 p 2 „ , 

Since the total number of tiles grows like — - — , this central region contains one half of the total number of tiles. 

By comparison, the situation is completely different in the other regions Ad', with D' < D. In this case, when D 
goes to infinity, the fraction of tiles in such a region vanishes. Let us first explicitly compute the asymptotic behavior 
of this fraction and then discuss how to handle this vanishing character. 



FIG. 16: (left) around the central region Ad, there is a first crown of triangles /Cd-i (light gray), a second crown of kites K,d-2 
(medium gray), a third crown of kites K.d-3 (dark gray) and so forth; in this figure, D = 6. Right: the three kinds of kites 
which define a kite-shaped region K, D i (white): a big kite minus two small kites (light gray) plus a smaller kite (dark gray). 

Figure Hrjl Cleft) shows the geometry of the regions under consideration when D is finite. Around the central region 
Ad, there is a first crown triangular domains, then a second crown of kite-shaped regions, and as k increases, the 
regions Ad' he in concentric crowns around the central one. Each crown contains 2D kite-shaped (or triangle-shaped) 
regions. As figure 1161 (right) illustrates, such a kite can be seen as a greater rectangular kite^S minus two equal 
rectangular kites. These two latter kites have a non empty intersection, which is itself an even smaller kite. If the 
areas of these kites are respectively denoted by 8d> , 6d'+i & n d Sd>+2, then the area of ICd 1 is 

kd< = 5 D ' - 2<5d'+i + 5 D '+2- (A. 9) 

Now, for each of these 4 kites, its edges are two rays of the circular central region and two tangents to this central 
region, perpendicular to these two latter rays. Since the central region has diameter p, if the two rays make an angle 
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6* the area of the kite is 



6(0) = ^ tan I- 

For the domain K.d>, = tt(D — D')/D for the greater kite. Henceforth, 

~ 2 f D-D n 



>D' 



P 



tan 7r 



2D 



The same formula holds for 8d'+i and 6d>+2- Thus if we set k = D — D' , 



p" 



' an ( -^-j^ ) 2 t an 



(k- 1)tt 
2D 



tan 



{k - 2)tt 
2D 



Fixing the ratio x = D' / D, and approximating the above as a second derivative, we finally get 

p 2 1 

*D> = k E (xD) ^ t -g^j tan (vr(l - ac)/2) [l + tan 2 (tt(1 - z)/2)] . 
Thus the number of tiles in Ae(xD) is 

N E (xD) = a E (xD)4>D(x) = 2Dke(xD)4>d{x). 
For later reference we also calculate the average area of a tile in Ae( x d), 

Ek, 



a(x) 



sin 2 cos (27ra:) - 1 + 2tt 2 x 2 



<i,j<E(xD) Sm D 



(TTX — Sm7TX 



Finally, the fraction of tiles in the region Ae(xD) is 

1 



nE(xD) 



Ne(xD) _ N E (xD) 



N T D 2 p 2 /2 2D 



tan(7r(l-a;)/2) [l + tan 2 (tt(1 - x)/2)] [ nx — sin Tix] 



(A.10) 



(A.ll) 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



(A.16) 



when D goes to infinity, since the number of tiles in the whole tiling is Nt = D 2 p 2 /2. 

As expected, these fractions vanish but they can be added to get a non vanishing number of intersections. More 
precisely, let us fix a real number < g < 1 and let us compute the fraction of tiles which lie in regions where at least 
E(gD) families of lines intersect: 



D-l 
l=E(gD) 



1 1 

2 + 2 



dxtan (irx/2) [l + tan 2 (ttx/2)] [tt(1 - x) 



sin7ra] = i + i/(l-0). (A.17) 



The above integral is a continuous function f(g) for g e [0, 1]. It fulfills the required condition /(0) = 1 (which means 
that the whole tiling contains a fraction 1 of tiles!). Moreover, f(g) < 1 when g > 0. Thus for any h < 1, there exists 
a real number g > such that a fraction h of the tiles lie in regions where at least E(gD) families of lines meet. Since 
g > 0, when D goes to infinity, the number of families of lines also goes to infinity in such regions. 

Moreover, appendix iBl demonstrates that the local tilings in such regions are true large codimension ones. Hence 
everywhere in the above region, the local entropy tends to the free boundary entropy a^. A fraction h of the tiles lie 
in regions where the local entropy per tile tends to when D — > oo. We gelpii the expected relation (|A.2|I and we 
conclude that free and fixed boundary entropies coincide. 

These arguments can also be extrapolated to finite p tilings provided D is large. Indeed, if p is finite, a domain ICd> 
can be very small and can contain very few tiles (and even no tile at all). Nevertheless in the previous demonstration, 
regions Ae( x d) with xq < x < xo + Sx can be put together into larger regions, to which all the previous arguments 
can be applied. 

Finally we remark that the tiling 7o which was not assumed a priori to maximize the entropy functional indeed 
does so at the large D limit. Therefore it should be close to generic and should give a good idea of the macroscopic 
structure of large D generic tilings, in particular as far as the scaling laws of section TA 31 are concerned. 
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2. Non-diagonal tilings 

For large p non-diagonal tilings with side lengths ki — Xip, where a < Xi < (3, we use a particular tiling T ' which 
is a variation of Tq. Its overall definition is the same except that the line spacing 1^ depends on the de Bruijn family 
-Fj. We follow the main steps of the prior demonstration, in particular the calculation of the areas K£>< or a pi. 

We choose h = p/(xip — 1) ~ l/xi so that the width U{ki — 1) = p of the family Fi is independent of i in the grid 
representation ffigure H31 left). Therefore the areas apr remain unchanged as compared to the previous section. 

On the other hand, the density of intersections in each domain Krv will vary because of varying line spacing. In 
particular the area ao of the unit cell now depends on i and j: 

= i hl i i = — r^in- ( A - 18 ) 

| sm Uij | XiXj | sin f/y | 

The number of intersections (|A.6J| becomes 

Qd>=kd>1 XiXjSin^-j^- I < /3 2 Qd>- (A.19) 

Henceforth, the new numbers of tiles in the crown Ad' and in the whole tiling satisfies Nd' < (3 2 Nrj> and Nt > a 2 Nr- 
If we set again D' — E(xD), the new fraction of tiles in Ad' satisfies 



N E (xD) . 1 , . 

n E (xD) = — r T < -o n E(xD)- A.20 

Nt a 



As a consequence, the outer crowns with small effective dimensionality D' , which had a vanishingly small contribution 
to the total number of tiles in the previous section, will again be negligible in the present case. 

To finish the proof, we must check that the remainder of the tiling has a local entropy equal to . The proof is 
explicited in appendix iBl in the diagonal case but can easily be adapted to the present case, leading to 

•--T- (A-21) 



(n% ~ P 2 7T 2 ' 

where the notations are defined in appendix [5] Even though the lower bound on tile fractions now depends on a and 
[3 because we are dealing with a nondiagonal case, it remains finite. 



3. Scaling law for the effective dimensionality 

We derive the effective dimensionality D' = D e ff(r) in the tiling 7o of side length p at the large D limit. Here 
< r < R is the distance from the tiling center O and R ~ pD/ir is the radius of the tiling. Assume that in a large D 
tiling, a region of effective dimensionality -D e // is an annulus Ad cJS of center 0, the tiling center, and of radii r and 
r + Sr. Hence D e ff(r + Sr) = D e ff(r) — 1 (the effective dimensionality decreases with increasing r), that is to say 

Sr^De ff {r) = -1. (A.22) 

Now we need to compute Sr. Following appendix I A II we write x = D e ff/D. If we knew the area s(x) of Ad b}} in 
the tiling metric, then we could extract Sr from s{x) = 2-nrSr, in order to write 

S (x)^-D eff (r) = -2wr. (A.23) 

But we can easily extract s(x) from equations i|A.14fl and l|A.15(l for the number and area of tiles in Ae( x d) to find 

Dp 2 

s(x) = — ^-tan(?r(l -x)/2) [l + tan 2 (tt(1 - x)/2)] [cos(2 nx) — 1 + 27r 2 x 2 ] . (A. 24) 



It is convenient to introduce dimensionless variables which remove the divergences for large D and p. Specifically, we 
introduce f = r/R and s(x) = s(x)/(Dp 2 ), and we recall that D e ff(r)/D = x and R = pD/ir to obtain 
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which is completely equivalent to eq. (|A.22I) . 

We can solve eq. (|A.25|) by direct integration to obtain the radius r corresponding to a given effective dimension x, 

r x i wx 

f 2 (x) - 1 = -n / s(x')dx' = — (-14 + 2cos(7rx) +7rxcsc 2 ( — )(irx + 2sin7ra)) (A.26) 
Jo 16 2 

or we can invert the solution and obtain (after returning to unsealed variables) 

D eff (r) = D 7 (r/R) (A.27) 

with 7(r) the inverse of r(x). Although the solution in eq. I|A.26[) cannot be inverted in closed form, we can expand 
it for small x to find values of r close to the boundary R. Defining the small quantity e = 1 — r/R, we have 

2 2 

+ 0(x 4 ), (A.28) 
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then by reversion of series obtain 
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7 (r/i?) = x « — sfe. (A.29) 



Apparently, the effective dimensionality varies rapidly near the boundary. 

We also use the notion of effective codimension c e // = D e ff — 2. At large D, one has 

Jy^~2 ~ ~^d~ = 7(r/i?) ' (A - 30) 

All the quantities expressed as functions of x are also functions of r/R. For example, the vertex density 

12 2 1 

d(x) = — — ~ — = — -— — ~ ; = d(r/R) (A.31) 

V ' <t(x) nx nj(r/R) ^6(1 -r/R) K ' 1 K 1 

at small x, i.e. when r close to R. This scaling compares well with numerical data in section Hvl 



APPENDIX B: REGULARITY OF LARGE CODIMENSION FIXED BOUNDARY TILINGS 

We check here that in each region Ad 1 of the diagonal tiling Tq the local entropy per tile tends to the free boundary 
one (Too as D — > oo. As discussed in paper li, a sufficient condition is that tile fractions should be "bounded" on Ad 1 - 
This means that all local tile fractions have the same order of magnitude as the strain-free ones n*- . In other 
words, there exist a finite constant a > such that 



Mj > a 



(B.l) 



for all i and j. We want to calculate the tile fractions n'^ in Ad 1 and compare them to the strain- free ones. Recall 
that in a D — > 2 tiling, 



n*. = C sin 



D 



(B.2) 



where the constant C comes from the normalization relation n*j = 1 . 

As in the previous appendix, we set D' = E(xD) where D is large and x > is finite. The main difficulty in 
this appendix comes from the fact that in the region Aw under consideration, not all tile species occur. Without 
loss of generality, we assume that only the tiles with indices i,j = X,...,D' exist. These tiles appear in the region 
Ae(xD) with the fractions i|B.2fl . but with a different normalization constant because only some of them appear. The 
normalization relation now reads 



E < 3 = l - 

l<i<j<D' 

At the large D limit, this sum can be replaced by an integral and we get the tile fractions in Ae( x d)'- 



(B.3) 



1 



(D 1 ) 2 ttx — sin(7rcc) 



sin 



J - * 



D' 



(B.4) 
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Now we check that this tiling has bounded local fractions of tiles when it is considered as a D' — > 2 tiling. We need 



to compare the above n'- to the corresponding quantities in an unstrained D' — > 2 tilings, in other words to 



ij 



J 



D' 



By convexity of the sin function on the interval [0, 7r], we have: sm\xir(j — i)/D'\ > xsm\n(j — i)/D'\, and 

6 



> 



> 



(n')*. 7r ttx — sin(7ra;) n 2 ' 



(B.5) 



(B.6) 



which achieves the proof: a = 6 /7r 2 in condition l)B.ljl . From the results of paper li, we conclude that in such a region 
Ad 1 , the local entropy per tile equals the free boundary value when D — > oo. 
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TABLE I: Exact (except B\\/B\q) data for fixed boundary tilings, p = 1 and D = 1, 2, . . . , 10. 



J) 


I 


2 


3 


4 


5 


hd\v — i) 


i 


i 

_L 


o 
z 


c 
o 




or> {p = 1) 







0.231 


0.347 


0.413 


Bd+i/ ' Bd 


1 


2 


4 


7.75 


14.65 


D 


6 


7 


8 


9 


10 


Bn(p=l) 


908 


24,698 


1,232,944 


112,018,190 


18,410,581,880 


od{p = 1) 


0.454 


0.482 


0.501 


0.515 


0.525 


Bd+i/Bb 


27.20 


49.92 


90.85 


164.35 


295.97 ±0.04 



TABLE II: Exact data for fixed boundary tilings, p = 2 and D = 1,2, ... ,6. 



D 


1 


2 


3 


4 


5 


6 


B D (p = 


2) 


1 


1 


20 


5,383 


16,832,230 


327,674,047,931 


&d(p = 


2) 







0.250 


0.358 


0.416 


0.451 



TABLE III: Convergence of D — 9 path count data for increasing run length. Measured numerical errors on the latter quantity 
(Num. err.) are always smaller than estimated error bars (Estim. err.). The latter are calculated using numerically measured 
standard deviations and conjectured auto-correlation times (section flV II [24]). 



Nmc 


10 3 


10 4 


10 5 


10 6 


10 7 


exact 


h 


159.264 


164.176 


164.411 


164.371 


164.344 


164.35 


logP 9 /9 


0.563396 


0.566771 


0.566930 


0.566903 


0.566884 


0.566891 


Num. err. 


3.5 10~ 3 


1.2 10~ 4 


3.9 10~ 5 


1.2 10~ 5 


6 10~ 6 




Estim. err. 


8.5 10" 3 


2.7 10" 3 


8.5 10 -4 


2.7 10" 4 


8.5 10" 5 
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TABLE IV: The first values of the limiting distribution of coordination numbers Z, obtained both in the mean-field approxi- 
mation and numerically, by Monte Carlo simulations. 



z 


3 


4 


5 


6 


7 


Mean-field 


0.33 


0.41 


0.20 


0.05 


0.009 


Numerical 


0.31 


0.43 


0.21 


0.04 


0.005 



TABLE V: Path count statistics for strained tilings. 



D 


p = 1, 1, ... 


p = 1,1, ...,2,2,... 


p= 1,1, ...,3,3,... 


p= 1,2,1,2,... 


10 


0.569 


0.572 


0.572 


0.573 


20 


0.573 


0.572 


0.571 


0.574 


30 


0.572 


0.571 


0.570 


0.572 


40 


0.572 


0.570 


0.570 


0.571 


50 


0.571 


0.571 


0.570 


0.570 



